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ABSTRACT 


A  unified  efficient  formulation  for  the  static,  dynamic,  and 
stability  reanalysis  of  locally  modified  structures  is  presented. 

The  reanalysis  problem  is  to  find  the  structural  responses  when  some 
of  the  element  properties  are  adjusted  as  a  result  of  design 
modifications.  The  reanalysis  problem  is  formulated  as  a  problem  of 
much  lower  order  than  the  original  system.  This  is  achieved  by 
utilizing  the  linearity  property  of  the  structure  using  the  pseudo¬ 
load  concept,  together  with  the  solution  of  the  original  system. 

The  modifications  to  the  structure  are  treated  as  displacement 
dependent  pseudo-loads  of  the  system.  By  expressing  the  modified 
system  response  as  linear  combinations  of  the  response  of  the 
original  system  and  a  term  depending  on  the  pseudo-load,  a  reduced 
set  of  response  equations  can  be  obtained.  In  the  static  and 
sinusoidal  steady  state  analyses,  this  leads  to  a  set  of  linear 
algebraic  equations.  For  free  vibration  and  stability  analyses, 
this  results  in  an  eigenvalue  problem.  General  modifications 
involve  multiple  parameters  so  that  numerical  solutions  are 
required.  Numerical  examples  are  included.  For  the  special  case  of 
a  single  parameter  variation,  a  simple  closed  form  formulation  of 
the  reanalysis  problem  is  provided.  Explicit  reanalysis  formulas 
are  presented  for  changes  in  the  cro3s-sectional  area  of  truss 
members,  the  stiffness  of  beams,  the  thickness  of  plane  stress 
elements,  and  lumped  springs  and  masses. 

INTRODUCTION 


Finite  element  analysis  is  an  emerging  technique  which  is  used 
in  the  design  of  complex  structural  systems.  From  spacecraft  to 
ground  vehicles,  from  the  largest  dam  to  a  small  welded  joint, 
almost  all  engineering  structures  have  been  analyzed  in  the  design 
stage  by  the  finite  element  method.  The  versatility  of  the  finite 
element  method  has  made  possible  the  analysis  of  complex  structural 
systems.  The  method  provides  the  capability  of  modeling  in  great 
detail  the  real  system  under  investigation. 

As  with  many  things  in  life,  this  increase  in  analytical 
capability  also  has  a  price  tag:  the  cost  of  setting  up  and 
executing  a  finite  element  analysis  and  the  loss  of  insight  into  the 
variations  in  the  structural  response  when  certain  design 
parameters  are  changed.  Both  of  these  obstacles  are  related. 

It  is  these  problems  that  jure  addressed  in  this  paper. 


S';.?  ' 

-need 


Ot 


*  va :  2 


cVa.T 


ft 


’ll 


and/,r 


;  ’tinl 


The  cost  problem  arises  because  in  the  application  of  finite  element 
analysis  to  design  many  configurations  are  investigated  before  the 
design  is  finalized.  These  investigations  are  a  natural 
consequence  of  the  trial  and  error  process  of  current  design 
procedures.  Also,  to  obtain  a  parametric  study  of  structural 
responses,  the  structure  must  be  analyzed  many  times  by  varying  one 
or  more  design  parameters.  Thus,  to  reduce  the  cost  of  finite 
element  analyses,  it  is  desirable  to  have  a  technique  which  will 
analyze  the  effects  of  system  modifications  efficiently  and 
economically,  especially  when  only  a  small  portion  of  the  original 
system  is  modified. 

Many  authors  have  treated  the  reanalysis  problem.  Reviews  of 
reanalysis  techniques  can  be  found  in  Wang,  Pala2zolo,  and  Pilkey 
(1)  and  in  Arora  (2).  In  particular,  the  static  reanalysis  problem 
has  been  investigated  by  Argyris,  Bronlund,  and  Roy  (3),  Kavlie  and 
Powell  (4),  and  Noor  and  Lowder  (S),  among  many  others  (1),  ( 2 )'. 

The  eigenvalue  or  free  vibration  reanalysis  problem  has  been 
studied  by  weissenburger  <6)  (7),  Stetson  and  Palma  (8),  Hallquist 
(9),  and  Hirai,  Yoshimura,  and  Takamura  (10).  The  forced  sinusoidal 
response  as  a  function  of  structural  parameters  has  been  studied  by 
Done  and  Hughes  (11)  and  verified  experimentally  by  Done,  Hughes, 
and  Webby  (12). 

The  purpose  of  this  paper  is  to  present  a  unified  formulation 
for  the  static,  free  vibration,  and  forced  sinusoidal  response 
analysis  of  a  locally  modified  structural  system.  The  structure  is 
assumed  to  be  modeled  by  finite  elements .  The  reanalysis  problem 
is  to  find  the  structural  responses  when  some  of  the  element  pro¬ 
perties  are  changed  as  a  result  of  design  modifications.  Because 
the  modifications  are  localized,  the  reanalysis  problem  can  be 
formulated  as  a  problem  of  much  lower  order  than  the  original  system. 

This  formulation  is  achieved  by  exploiting  the  linearity  property  of 
the  structure  using  the  pseudo-load  concept,  together  with  the 
solution  of  the  original  system.  The  modifications  to  the  structure 
are  treated  as  displacement  dependent  pseudo-loads  of  the  system. 

Depending  upon  the  type  of  analysis,  these  displacements  (fl)  may 

be  parts  of  the  static  displacements  or  the  natural  vibration  mode 

shapes  of  the  system.  By  expressing  the  modified  system  response  as 

linear  combinations  of  the  respones  of  the  original  system  and  a  term 

depending  on  the  pseudo  load  (and  hence  on  (fl),  a  reduced  set  of 

equations  with  fl’  as  the  unknown  can  be  obtained.  In  the  static 

and  sinusoidal  steady  state  analyses,  this  procedure  leads  to  a  set  of  linear 

algebraic  equations  in  fl'.  For  free  vibration  analyses,  the 

procedure  results  in  an  eigenvalue  problem  in  the  form  of  an  r 

order  determinant,  where  r  is  the  number  of  degrees  of  freedom 

affected  by  the  modification.  For  this  case,  if  all  the  modes  are 

used  in  the  formulation,  the  resulting  equation  is  the  exact 

characteristic  equation  of  the  modified  system,  otherwise,  an 

approximate  formulation  is  obtained.  Methods  of  improving  solution 

accuracy  when  only  a  partial  set  of  modes  is  available  will  be 

discussed.  General  modifications  involve  multiple  parameters  so 

that  numerical  solutions  are  required.  For  the  special  case  of  a 

single  parameter  variation,  a  simple  closed  form  formulation  of  the 

reanalysis  problem  can  be  derived.  A  few  explicit  reanalysis 

formulas  are  presented  in  the  following  sections. 


GENERAL  REANALYSIS  PROBLEM;  PSEUDO-LOAD  FORMULATION 

Consider  a  general  linear  structural  system  which  can  be 
described  by  the  following  equation 

D  U  -  f  (1) 

where  D  is  the  generalized  stiffness  matrix.  For  a  static  problem, 

D  -  K 

For  a  steady  state  sinusoidal  problem 
D  *  K  -  cj2M 

Suppose  now  the  system  is  modified  so  that 

D'  -  D  +  AD  (2) 

Substitute  (2)  into  (1)  and  rearrange 


DU'  -  L  -  AD  U' 


(3) 


where  u'  is  the  solution  of  the  modified  system. 
Fremultiply  (2)  by  2.  the  inverse  of  0, 


il’  -  2  1  -  2  AD  11'  (4) 

For  local  modifications,  bp  will  be  a  highly  sparse  matrix; 
that  is,  Ap  contains  only  a  small  number  of  nonzero  columns  and 
rows .  Let 


AD 


a  submatrix  of  AD  which  contains  only  the  nonzero 
columns  of  AD 


Furthermore,  let  the  nonzero  column  of  AD  be  the  j  ,  j 
j  columns,  and  let  * 

u 

*1 


u' 


u. 

32 


u. 

3, 


9  • 


ug  ■  the  column  vector  formed  by  deleting  from  u* 

Similar  definitions  apply  to  A  and  ug.  With  these  definitions, 

(1)  can  be  simplified  to 


/N  A 

li'  ■  U  “  2  AD  U 


Now,  define 


/N 

Xi  .  Z  ADf 


X  •  Z  AD 


,th  ✓N 

where  AD^  “  i  column  of  .AD,  Equation  (_5 )  can  be  written  as 


u'  -  u  -  Y  u1 


Rearrange  the  equations  in  (2)  to  get  the  following  partitioned 
set  of  equations. 


From  the  first  r  equations  of  (fl). 


A.  A  A  A, 

u  ’  -  u  -  Y  u ' 


A  A  -1  A 

U*  -  ( I  +  Y)  U 


and  <J’g  can  be  computed  from  the  lower  portion  of  equation  (fl). 


A.  A  A  A, 

u '  -  u  -  Y  U  ' 
-S  -9  -s  - 


Note  that  equations  (3)  and  (IQ)  provide  the  solution 
of  the  modified  problem.  The  computational  efficiency  obtained  in 
solving  the  modified  system  formulation  stems  from  the  fact  that  equation  (3 ) , 


f 


which  is  a  system  of  r  equations,  is  to  be  solved  instead  of  solving  n 
equation*,  which  would  be  the  case  if  the  modified  system  were  to  be 
solved  directly.  The  additional  information,  i.e.,  matrix  Y,  needed 
in  this  formulation  can  be  obtained  by  solving  the  original  system  with 
r  different  loadings .  This  solution  is  accomplished  by  forward-backward 
substitution  once  the  original  general  stiffnes  matrix  is  factored 
into  upper  (H)  and  lower  (&)  triangular  forms. 

Note  that  in  the  abo^e  derivation  we  used  only  the  linearity 
property  of  the  system.  Alternatively,  it  can  be  shown  that 
equations  (3)  and  (lfl)  can  be  derived  directly  from  the  Woodbury  formula  of 
numerical  analysis,  which  is  also  known  as  the  matrix  inversion  lemma 
in  the  control  theory  literature. 

In  the  following  sections,  equations  ;S)  and  (IQ)  will  be  applied 
to  the  reanalysis  of  static  response,  steady  state  sinusoidal  response, 
and  free  vibration  problems. 

STATIC  REANALYSIS 

Por  a  static  problem 


D  -  K 


(11) 


and  Z 


is  the  system  flexibility. 


Since 


AD  -  AK 


(12) 


equation  (£a)  becomes 


-1  /N 

-  K  AK.  (13) 

where  is  the  ith  nonzero  column  of  AK.  Assuming  the 
original  system  is  solved,  then  u  is  known,  and  &  can  be  written 


K  -  L  U 


(14) 


In  terms  of  upper  and  lower  triangular  forms,  equation  (l)  becomes 


U  z  -  f 


(15) 


L  U 


z 


(16) 


The  reanalyais  procedure  is 


1.  Formulate  the  modification  in  the  stiffness  matrix  as 


n 

AK 

x-1  e. 


(17) 


where  AX  •  the  contribution  in  the  global  stiffness  matrix  due 
e  th 

i  to  the  modification  of  the  i  element  stiffness 

characteristics . 


2.  compute  matrix  Y  by  the  following  algorithm.  The  i 
column  is  computed  as 


th 


U  2, 


✓N 

(18) 

z . 

- 1 

(19) 

where  A **  i  nonzero  column  of  AK.  In  equation  (18),  . 

2 .  is  computed  by  backward  substitution;  and  from  (12) 

X.  is  computed  by  forward  substitution .  Carry  out  these  procedures 
in  (Ifl).  (12)  for  i  -  1,  2,  ...  r. 

.  .  .  .  .  .  . .  .  th  .  th  th 

3.  Form  matrix  X  by  selecting  the  j  ,  j  »  •  •  •  jr  row 

of  X  •  '  1 

4.  Compute  fl’  from 

A  A 

W  u’  -  u  (20) 

where 

W  -  I  +  Y  (21) 

Equation  (2S2)  can  be  solved  by  decomposition; 


W  -  L 

U 

(22a) 

-w 

-w 

■ 

IN> 

A 

-  u 

(22b) 

A. 

A 

(22c) 

L  u’  - 
-w  - 

z 

5.  Compute  ji'  from  equation  (2)- 

This  completes  the  solution  of  the  modified  system.  Note  that  in 
the  final  step,  fl'  will  be  computed  again.  The  results  can  be  used 
to  check  against  2’  obtained  from  the  solution  of  equation  ( 20 ) . 

Note  also  that  in  this  computational  procedure,  the  matrix  partition  indicated 
in  equation  (fl)  does  not  need  to  be  carried  out  explicitly. 

For  general  modifications  involving  many  degrees  of  freedom  (dof), 
the  above  procedure  can  be  programmed  as  a  post-processor  of  any 
general  finite  element  code.  For  simple  modifications  involving 
few  degrees  of  freedom,  the  reanaly3is  equations  can  be  formulated 
manually  from  finite  element  solution  outputs.  It  is  also  possible 
to  develop  a  "closed  form"  formula  for  simple  local  modifications. 

To  illustrate  the  development  of  an  analytical  modification 
expression,  consider  the  case  of  adding  spring  Ak  between  dof 
and  dof  j  in  a  general  finite  element  system. 

Step  1:  For  this  case,  the  first  step  of  the  procedure  is 
readily  accomplished  since  it  is  known  that 


AK 


0 


0 


Ak 

-Ak 


0 


j 


1 


-Ak 

Ak 


^  2 


0 


j2 


(23) 


That  is,  AK  has  only  four  nonzero  elements.  It  follows  that 


r  o  'i 


A 


AK 


0 

+Ak 

0 

«  :  ► 

o 

-Ak 

0 


is  the  first  nonzero  column  of  AK,  and 


0 


Step  2 1  Compute  matrix  Y  using 


&  11  -  AKX  -  Ak  £ 


K  i2  •  4K2  -  -  ik  £ 


(25a) 

(25b) 


For  this  special  case,  it  is  sufficient  to  solve 


K  Y  -  P 


(26) 


This  can  be  accomplished  using  the  procedure  of  equations  (iS)  and  (13). 
after  obtaining  Y^, 


and 


Y,  -  Ak 

Y, 

-1 

-1 

Y„  -  -Ak 

Y_ 

-2 

-2 

i  -  : 

Y2]  *  Ak 

r  y  -  y 

~1  -2 

] 

27) 

3:  Form  which 

is  the  matrix  formed  by  the  j  th  and 

row  of  Y.  Hence 

X 

V 

V 

A 

Y  *  *k 

V 

•V 

(28) 

Step  4:  To  find  ; 

&  \  first 

form  M  as 

1  +  b 

-b 

W  -  I  +  Y  - 

(29) 

a 

1  -  a 

where 


a 


Ak  Y 


b 


Ak  T 


Hence 


which  can  be  expressed  explicitly  as 


-  C 


(1-a) u.  +  bu . 

31  32 


-au.  +  (1+b) u . 

3,  3- 


where  C  -  1  +  b  -  a.  From  (IQ)  and  (il) 


(30) 


(31) 


(1-Ak  r.  )  u  +  Ak  7  u. 
V  J1  -3X1  32 


c 


1  + 


1 


)  Ak 


(32) 


f 


Note  that  for  this  case,  we  have  obtained  (and  hence 
y)  explicitly  as  a  function  of  Ak,  the  varying  parameter. 

As  a  consequence,  we  have  solved  a  family  of  problems,  as  well  as 
a  particular  reanalysis  problem,  since  for  any  value  of  Ak,  the  response 
of  the  modified  system  can  be  computed  by  simple  calculations  using 
equations  (22)  and  (2). 

As  a  numerical  example,  consider  the  simple  truss  shown  in  Pig.  l. 
Suppose  now  that  the  cross-sectional  area  of  member  4  is  changed  to 
A^  -  (  V2/2 )( 1+ar ) .  we  want  to  find  the  response  of  the  modified  system. 
The  original  system  stiffness  is 


K  -  0.125(106) 


5  1 

1  5 

0  0 

0  -4 


0  0 

0  -4 

5  -1 

-1  5 


and  the  decomposed  matrices  are 


K  -  L  U 


P,  =  1000  lbs 


L2  *  L4  = 

*1  =  *3  = 

A2  “  A4  " 
E  =  1.0  x 


(4448. 2N) 

Lg  =  20  in.  (50.8  cm) 

20 in.  (71.84  cm) 

Ag  =  1.0  in2  (6.45  cm2) 

ST/  2  in2  (4.56  cm2) 

7  2  -> 

10  lb. in  (68.45  GN/m") 


where 


Fig.  1  Truss  used  as  an  example  problem 


L  -  0.125(10®) 


V  5 
1V5 
0 


0  0  0 

V24/5  0  0 

0  V5  0 


0 


-V3/10  -1/V5  V22/5 


u 


0. 125(106) 


The  loading  vector  is 


The  displacements  o£  the  original  system  are  calculated  to  be 


r 


u  -  (2/11)10 


-3 


For  the  modified  system,  it  can  be  shown  that 


AK 


0 

0 

0 

0 


0  0  0 

0  0  0 

0  e  -e 

0  -e  e 


where  e 


0.125( 


10 


6 


)0t. 


Thus,  the  matrix  Y  can  be  computed  to  be 


1  -1 

-5  5 

Y  *(«/ll)  1  -1 

-6  6 


Wow,  to  use  equation  (12).  we  have  ^  •  3,  j2  ■  4  and  calculate 


(10/11) (10‘3) 


Uj  -  u4  -  (50/11) (10-3) 


Y.  -  1 
^l1 


V  ■  -6 

Ak  =  a/ll 
From  equation  (12) 


C  -  1  +  (Yi  i  -  Yj  .  )Ak  **  1  +  7a/ll 
J1  2 

U'  -  (1/C)[(1  +  6a/ll)u3  +  ( a/ll) ] 
-  (1/C)  [6au3/ll  +  (1  +  at/ll)u4] 

and 


-  ux  -  a(ii-  -  Jl^/ll 

-  5a(-u3  +  11^ ) / 1 1 

Hence,  the  modified  system  solutions  axe  found  in  terms  of  the 
original  system  solution  and  the  design  parameter. 

Note  that  equations  (13)  offer  an  opportunity  to  select 
design  parameter  a  for  prescribed  static  responses.  For  example, 
from  equations  (33a  and  b),  we  can  solve  for  a  as 


a  -  11(U3  -  Uj)/(7u^  -  6u3  -  U4) 

Suppose  we  want  the  modified  system  to  have  a  smaller  deflection 
at  dof  3.  Say  we  want  u'3  “  0.9  u3*  Then  from  equation  (14), 
the  required  a  is 


( 33a) 

(33b) 

(33c) 

(33d) 

( 33e) 

(34) 


a  -  0.234 


Thus,  using  the  modification  theory,  we  can  perform  a  parametric 
study  for  simple  design  modifications,  as  well  as  design  the  system 
for  prescribed  responses  by  properly  choosing  the  design  parameter 
value . 

Although  only  one  special  case  of  adding  a  spring  between  two 
dof  has  been  treated  analytically  here,  many  other  special  cases 
can  also  be  derived  in  a  similar  manner. 

STEADY  STATE  RESPONSE  REANALYSIS 

For  a  steady  state  analysis,  let  the  forcing  function  and  the 
responseble  sinusoidal  functions  of  time,  i.e., 

f  ■  f  sin  ut  and  u  *  u  sin  ut 


where  oris  the  forcing  function  frequency, 
systems.  The  general  stiffness  matrix  is 


This  applies  for  undamped 


K  -  w  M 


(35) 


The  modification  will  involve  changes  in  both  the  mass  and 
stiffness  matrices.  Thus 


AD 


AK  -  (j  AM 


(36) 


If  this  frequency  response  problem  is  to  be  treated  directly,  the 
modified  system  can  be  solved  using  the  pseudoload  formulation  with 
equations  (3S)  and  (36)  for  D  and  AD.  If  the  modal 
superposition  method  is  used  for  the  analysis,  then  the  receptance 
matrix  Z  cam  be  represented  using  a  modal  expansion 


(37) 


P 

Z 

4=1 


4  $ 


^  (ui  ~  "*> 


where  w.  is  the  frequency,  <t>  .  is  the  mode  shape,  and  G  is  the 
=  th  •  * 

generalized  mass  of  the  l  mode.  If  P  *  n,  equation  (22)  is 

exact.  If  p  <  n,  which  is  usually  the  case  with  the  normal  sinusoidal 
mode  method,  equation  (22)  gives  approximate  results.  The 
results  can  be  improved  by  isolating  the  static  response  in  Z.  it  can 
be  shown  that  equation  (22)  can  be  written  as 


n 


(38) 


z.  I  *■  * 

-  ^7^) 


z 

1=1 


R-l 

Z 
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4tf/U2Y 

v?w  / 


<w 
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Th#  integer  R  is  a  higher  mode  factor.  Normally,  the  higher  the 
value  of  R,  the  better  the  accuracy  and  lower  the  computational 
efficiency.  It  can  be  shown  that 


-1 


n 

1=1 


K-1 (M  K'1)3 


l  ^  ^ 

1=1  V^2)a+i 


(29) 

(40) 


Rewrite  (.18)  as 


-  E 

i  =1  V*?  \  ^  ' 


Using  (IQ)  and  (41), 


R-l 

2  =  Zn 
—  a=0 
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W)a 


Ll=i  (u£)a+ 
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R-l 

zn 

a=0 


(a)2)  3  (K-1)  (M  K_1)a 


\  /  a>2  \ 

1=1  «  Vo,*  / 


(42) 


The  first  series  on  the  right  hand  side  of  (j Z)  can  be  found  from 
the  static  solution  and  is  independent  of  the  number  of  modes 
available.  Bence,  the  improved  receptance  formula  when  only  p 
modes  are  available  is 


Z 


R-l 
a 


ZQ  (K_1) (M  K_1)d 


(43) 


•  To  find  u  and  2,  proceed  in  the  following  manner. 


u  ■  Z  f 


R-l 


E  u 

a*0  —a 


P 

E 

1*1 


(44 


where 


“a  * 


(45) 


Z  -  K  1(M  K_1)a 
a 

Compute  u  from 

K  u  -  f  (46a) 

-  -o 

and  calculate  u  using 
a 

K  u  -Mu.  for  a  -  1  to  R  -  1  (46b) 

- a  -  - ( a-1) 

This  procedure  can  also  be  applied  to  find  2 .  Once  1  is 

obtained,  one  can  proceed  as  in  the  static  reanalysis  problem  to 

find  the  response  of  the  modified  system  using  equations  (2)  and  (IQ). 

FREE  VIBRATION  REANALYSIS 

If  the  external  forcing  function  is  zero,  then  u  -  0  and  equation 
(2)  has  a  nontrivial  solution  if  and  only  if 


d«t  (I  +  $)  -  0 


(47) 


This  is  the  frequency  equation  for  free  vibration  of  the  modified 
system.  For  the  nodal  formulation,  the  matrix  Y  can  be  found  using 
the  formulation  of  the  previous  section.  For  free  vibrations,  equation 
(2)  becomes 


By  using  the  calculated  modified  system  natural  frequency  u.  in  X 
of  Eq.  (AS),  the  system  mode  shape  can  be  computed  from  (ifll  by 
properly  assigning  a  unity  value  to  one  of  the  components  of  u' .  Any 
component  of  y'  will  do  as  long  as  it  is  not  a  node  point  of  that 
mode. 

In  the  following,  we  will  develop  several  special  cases  of 
frequency  equations  of  modified  systems,  we  will  use  the  modal 
superposition  approach  assuming  all  modes  are  available  and  then 
indicate  how  to  improve  the  results  when  only  a  partial  set  of 
modes  are  to  be  used. 

Consider  the  special  case  of  adding  mass  Am  an  dof  j.  In 
general,  a  point  mass  in  a  three-dimensional  space  has  3  dof.  Here  it 
will  be  assumed  that  the  mass  is  constrained  for  one-dimensional  motion  only. 
Neglect  the  contribution  to  the  stiffness  due  to  this  added  mass 
(or  mass  removal,  if  Am  is  negative).  Then 


From  (56)  and  (47),  the  frequency  equation  for  adding  Am  at  dof 
jl  is  ~ 


uTAm 
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1=  1 


0  jfc  - 
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Let  w  -  u‘.  be  the  solution  of  equation  (52).  The  mode  shape  of  the 
modified  system  corresponding  to  natural  frequency  can  be 

‘computed  from 


♦  ■  *=  Y  (0)  =  CJ.‘ ) 

-l  —  x 


where  X  is  given  by  equation  (52).  In  the  above  derivation,  if 
the  approximate  Z  defined  by  equation  (42)  is  used  in 
equation  (52).  then  an  approximate  frequency  equation  is 
obtained . 

* 

Consider  next  the  case  of  adding  spring  Ak  between  dof  j  and 
dof  j  For  this  situation,  AD  equals  Ak  of  equation 
( 23 )  so  that 

/\ 

ADX  ■  Ak  E 
AD2  -  -Ak  E 

with  £  given  by  equation  (24£).  Thus, 
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From  equations  (££)  and 
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(64a) 


(64b) 


( 64c ) 
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Substitution  of  equations  (54)  into  (5^)  leads  to  the  frequency 
equation 


n 

1  +  Ak  Z 
)l-l 


“V  '  V* 


0 


(65) 


Equation  (65)  is  the  frequency  equation  after  adding  spring  bk  to 
the  system.  To  find  the  mode  shape  of  the  modified  system,  choose 
flL'  ■  1.0,  then 


1  +  Y 


V 


V 


(66) 


k 


Having  determined 


the  mode  shape  can  be  computed  using  equation  (2). 

Consider  now  the  special  case  of  adding  a  spring  between  dof 
j  and  a  fixed  base.  For  this  case,  <p .  -  0  and  the  frequency 

1  ^2  *■ 

equation  reduces  to 
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2.-1 
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Ak 


(67) 


Furthermore,  let  Ok  -  <*>,  that  is,  introduce  a  rigid  constraint  at 
dof  j  .  The  frequency  equation  becomes 
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Z 

2-1 
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G  y?  - 


(68) 


As  a  numerical  example,  consider  the  truss  of  Fig.  l  with  unit  mass 
at  nodes  A  and  B.  Then  for  free  vibration,  the  original  system  mass 
matrix  is 


CM] 
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0  0 
1  0 
0  1 
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(69) 


The  corresponding  stiffness  matix  is 


[K]  -  a 
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a  -  0 . 125 ( 106 ) 


The  natural  frequencies  are 
u*  -  0.7639a 


w. 


4.764a 


2 

U) 

3 


5.236a 


9.236a 


with  the  corresponding  mode  shape  vectors 
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'  0.6882 

'  0.1625" 
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0.1625 
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0.6882 

1  i3 
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It 

-ef 
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,  4 

(72 


Note  that  the  mode  shapes  have  been  chosen  to  yield  the  unit 
generalized  mass,  or  -  l.o  for  t  -  1  to  4. 

Consider  the  modification  of  doubling  the  cross-sectional  area 

2 

of  member  3,  that  is,  let  AA3  ”  1-0  in  .  Since  member  3  acts  like  a 

spring  between  dof  2  and  4,  we  can  use  equation  (&£)  as  the  frequency 

_  2  2  _  2  2 

equation  for  the  modified  system.  Define  acj  -  w  ,  au,  -  .  Use 

Ak  -  AAjE/Lj,  *  2,  j2  -  4.  Equation  <fi5)  becomes 


(73) 
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Using  the  modal  data,  this  frequency  equation  becomes 


Q - 1056 _ 

4.764  -  u2 


_ LiM5_ 

9.236  -  J2 


-0.25 


(74) 


Equation  (24)  is  the  frequency  equation  of  the  modified  system.  It 
may  appear  odd  that  equation  (24)  has  only  two  roots  when  the  system 
has  4  dof.  Thi3  condition  results  because  in  mode  numbers  1  and  3  of 
the  original  system,  there  is  no  relative  motion  between  dof  2  and 
4,  and  hence  the  first  and  third  natural  frequencies  will  not  be 
altered  by  the  modification  of  the  cross  section  of  member  3.  The 
solutions  of  equation  (24)  are  the  new  frequencies  of  modes  2  and  4  of  the 
original  system.  These  solutions  are 

u2  -  4.917,  17.083 


which  are  the  same  as  in  the  direct  solution  of  the  modified 
system. 

It  is  of  interest  that  in  equation  (22),  we  can  specify  a  natural 
frequency  5  of  the  modified  system  and  compute  the  required  ■ 

This  procedure  provides  a  means  of  designing  for  a  specified  natural 
frequency. 

The  general  procedure  set  forth  here  can  be  used  to  derive 
frequency  equations  for  other  particular  cases,  such  as  a  change  in  the 
El  of  a  beam  element.  This  will  be  the  subject  of  future  presentations. 

CONCLUDING  REMARKS 

In  this  paper,  a  unified  formulation  is  presented  for  the  efficient 
solution  of  reanalysis  problems  for  locally  modified  systems. 

This  technique  is  applicable  to  both  static  and  dynamic  ( steady 
state  and  free  vibration)  reanalysis.  For  simple  reanalyses 
involving  one  parameter,  closed  form  formulas  are  derived  and 
procedures  of  using  these  are  outlined.  These  closed  form  formulas 
not  only  solve  the  desired  reanalysis  problem,  but  they  also  provide  a  ' 
means  of  performing  rapid  and  economical  parametric  studies  of 
structural  systems.  Furthermore,  for  one  parameter  problems,  the 
reanalysis  formulation  can  be  inverted  to  solve  the  problem  of 
selecting  a  parameter  value  which  results  in  a  given  system  response. 

For  modifications  involving  more  than  one  parameter,  the  formulation 
transforms  the  problem  into  one  equal  dimension  to  the  number  of 


degrees  of  freedom  affected  by  the  modification.  Thus,  it  would  appear 
that  this  formulation  can  be  programmed  on  a  small  de3k  top  (or  home) 
computer.  The  original  system  can  be  solved  on  a  large  computer  and 
then  a  small  computer  can  be  used  for  computing,  storing,  and 
transferring  the  required  data  for  reanalysis.  The  reanalysis  can  be 
performed  economically  and  rapidly,  if  the  tunr-around  time  as  well  as 
computing  time  is  considered,  for  many  local  modifications  during  the 
design  process.  Hopefully,  this  approach  can  also  furnish  more  insight 
insight  into  the  behavior  of  a  complex  structural  system  as  a  local 
parameter  varies. 
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